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Abstract 

We consider spinless electrons in two dimensions with the bare spectrum 
e (p) = \Px\ + \Py\- in momentum space, the interactions among electrons 
have a finite range qo, which is small compared to the Fermi momentum. 
A golden rule calculation of the electron lifetime indicates a breakdown of 
the Landau Fermi liquid in the model. At the one-loop level of perturbation 
theory, we show that the density wave and the superconducting instabilities 
cancel each other and there is no symmetry breaking. We solve the model via 
bosonization; the excitation spectrum is found to consist of gapless bosonic 
modes as in a one-dimensional Luttinger liquid. 



1 



Typeset using REVTpjX 



Recently, the question about validity of Landau's Fermi liquid (LFL) theory in two di- 
mensions has attracted much interest. This is because there are strong indications that 
in the normal state of the high temperature superconductors (HTSC's), some kind of a 
non-LFL state is realized. This idea is due to Anderson [[], who realized that there are 
striking analogies between the phenomenology of the cuprates and of one- dimensional met- 
als. The latter are known to be Fermi liquids of a very special type (Luttinger liquids): the 
Fermi surface (defined as a set of points, where the momentum distribution function has a 
singularity) still exists and has the same volume as for non-interacting electrons in agree- 
ment with Luttinger's theorem, but elementary excitations are completely different from 
those of a usual metal in D = 3: there are no quasiparticle excitations and the spectrum 
consists of collective modes only. However, in a D = 2 low density Fermi gas with short 
ranged interactions, the LFL phase was shown (within perturbation theory) to be stable 
0. It is believed, that LFL-behavior can break down in perturbation theory only in specific 
situations: assuming particular bandstructure (like in the nested Fermi liquid theory [[|) 
and/or singular long ranged interactions. The latter suggestion is due to Haldane j3j (see 
also ||), who proposed a generalization of bosonization as a tool to study non-LFL behavior 
in higher dimensions. Haldane's method was applied recently to study fermions interacting 
with long-ranged current-current interactions |J. 

Here we study a model with nesting. We consider spinless fermions in D = 2 described 
by the Hamiltonian 

# = 5>P C P C P + 7^ ^ ^(q) C P+q C P'-q C P' C P, 

p p,p',q 

where e p = vp {\p x \ + \p y \) is the bare electron spectrum, V^(q) is a Fourier transform of 
the electron-electron interaction and Q is the (two-dimensional) volume of the system [|7| . 
There are 4 branches a = 1,..,4 of the spectrum, one in each quadrant of the p plane 
(see Fig.l). The group velocity in branch a is v„ = de^/dp. Our dispersion corresponds 
to fermions moving in a planar quadratic lattice described by the hopping Hamiltonian 
H k in = Ei,R^RCiC i+R , where t-R,/v F = n5 mfi 5 nfi + 5 nfi ((-l) m -1) firm 2 + <5 m , ((-l) n -l)/vrn 2 . 
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We denoted R = (m,n); S m ,n is the Kronecker symbol. It is seen that the chosen band 
structure is quite unrealistic, since it requires that the overlaps of localized orbitals fall off 
only quadratically with distance. The motivation for studying the model Eq.(l) is twofold: 
the linearity of the spectrum allows for a solution of the model via bosonization; on the 
other hand, the Fermi surface of the model Eq.(l) exhibits the property of perfect nesting 
without being complicated by the occurence of van Hove singularities, as is the case for 
the nearest neighbor tight binding model. Thus it provides a convenient testing ground for 
nested Fermi liquid theory. 

Let us assume now that V(q) = VQ(qo — |q|), Q(x) is the Heaviside function. We take 
q <C A, where A is the length of the Fermi surface in one of the branches (see Fig. 1). 
For what follows, it is sufficient to consider potentials with finite range in momentum space 
and our choice is dictated by convenience only. In real space, our potential falls off at large 



distances q r > 1 as V(r) = (2^2 /n)U sin(g ?" — n/4)/(q r) 3 / 2 , where U = Vql/Ait stands 

for the contact interaction energy. 

We calculate the lifetime of an electron using Fermi's golden rule. There are three types 
of scattering possible: an electron p scattering to the state p + q in branch 1 can excite 
an electron-hole pair k and k — q in branch 3, branch 1 or branch 2 (4). We denote the 
corresponding processes by #4 and gpL- The Qfl processes do not make use of the 
peculiarities of the bandstructure and lead to the scattering rate 1/r w G 2 e 2 /E F of Fermi 
liquid type, where G = (1/2tt) 2 AV/v f is a dimensionless coupling constant. For the g 2 
processes, the scattering rate in the limit e, e p <C vp% becomes 



demonstrating the well-known linear dependence of the lifetime on energy. Consider now 
the <7 4 processes. Since the Fermi surface is a line and not a point, scatterings in this channel 
do not cancel as in the case of spinless fermions in D = 1, and there will be only partial 
cancellation of such scatterings due to Pauli principle. This happens if all momenta involved 
in the scattering process differ by q . In our case q <C A and cancellations are completely 
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^ -jr G 2 (e-e p )0(e-e p ), 



(2) 
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irrelevant. For e, e p <C VFqo, we find 

^=2vr|G 2 e ^( e -e p ), (3) 

where 5(x) is a Dirac delta-function. Eq.(3) is analogous to the scattering rate of electrons 
with spin in D = 1 in the g± channel ||. (In one dimension, this golden rule result indicates 
the spin-charge separation found in the exact solution ||. In fact, by Kramers- Kronig analy- 
sis one finds that the real part of the self-energy is 3?S(p, e) ~ 1/ (e — e p ). The 'quasiparticle' 
spectrum is given by the solutions of e — e p — E(p, e) = 0. There are two of them, corre- 
sponding to spin and charge excitations.) In the present case, the spin degree of freedom is 
replaced by the position within the branch and, as a result of g^-like scatterings, coherent 
propagation of electrons is again lost. A similar golden rule calculation was performed in 
the context of nested Fermi liquid theory of HTSC's ||. Unfortunately, the authors of Ref. 
H missed the contribution of g 4 -like processes to the scattering rate and, as a consequence, 
also the spin-charge separation. 

Due to the perfect nesting of the Fermi surface, the model Eq.(l) shows a strong tendency 
towards the charge density wave (CDW) formation. This state is characterized by a nonzero 
expectation value of the operator D p = c p+ qC p , describing a particle-hole pair with total 
momentum Q (where Q is the nesting vector of the Fermi surface, see Fig.l). It is important 
to note that there are in general two nesting vectors Qi and Q2 of the Fermi surface. We 
restrict ourselves to the case of a half filled band, for the moment. In that case Q2 = Qi + G 
where G is an inverse lattice vector and Q 2 = Qi = Q. This in turn means that an electron- 
hole pair D p can be scattered around the whole Fermi surface. Rewriting the interaction 
term of the Hamiltonian, the scattering in the CDW channel is found to be — V(p, p')-D p /-D p . 
Assuming that a charge density wave is formed, we can write (D p ) = d p (where (...) denotes 
the thermal expectation value) and the mean field equation for the order parameter d p reads 

^P = ^En P ,k)^-tanh^, (4) 

where E p = ^Je p + d p . Repeating the above analysis for the superconducting instability, we 
assume pairs of particles P^ = c p cL p to be scattered coherently around the Fermi surface. 
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The relevant matrix element is V(p, p') PL P p and the equation for the expectation value 
of the order parameter (Pj) = A p has the familiar BCS form 



with P p = ^/e p + A p . Note that the sign in front of the interaction term V(p, k) is different 
for the CDW and BCS instabilities. To determine the kind of instability which occurs (on 
the mean field level), one has to find all nontrivial solutions to Eqs.(4,5) and determine the 
corresponding T c values. The instability with the highest value of T c is dominant flO| . 

In order to decide whether the instabilities are real, we investigate two body scattering 
vertices in the corresponding channels and investigate their scaling in one loop approxima- 
tion. The mean field approximation corresponds to restricting ourselves to one particular 
one loop diagram. We will show that taking into account all diagrams of the given order 
completely changes the result and instead of a logarithmic flow, one observes that the scat- 
tering vertex stays constant suggesting a line of fixed points in our toy model, as is the case 
for D = 1 [II]. Let us start with the CDW channel. There are 5 one loop diagrams for 
scattering from Pp to PL: (a) the ladder diagram, (b) the crossed diagram, (c) + (d) vertex 
corrections and (e) the screening diagram. We set p and p' to lie on the Fermi surface and 
the energies of the incoming and outgoing electrons are 0. We work at finite temperatures. 
The diagram (a) is the one considered in the mean field theory. Its contribution is precisely 
cancelled by (b). In fact, the condition q <C A implies small momentum transfer scatterings 
and that is why p « p' rs k. For such momenta, an exact equality e_k+ P + P ' = — ek holds 
and, consequently, 

r (o) _ r (6) _ V- v rr /fck) - /(-Ck) 

1 p,p' - _i p,p' - 2^ V pM v Kp' • W 

The remaining diagrams (c-e) contribute 

r p c pf +e V) = v p>p/ £ ww /(ekW(ek+p, - p) , (7) 

k Ck — Ck+p'-p — w 

where W p p : ^ = — V p p i + Vp+ q ,k + Vp,k and we have introduced an infinitesimally small 
energy transfer u into the scattering process. Evaluating T p c ^f +e ^ (a?) for finite u and taking 



the limit u — > gives Tptf +e ' = 0. Thus the scattering vertex in the CDW channel does 
not flow in the one loop approximation. The analysis of the BCS channel is analogous. 
Summarizing, we do not find any sign of symmetry breaking neither in the CDW, nor in the 
BCS channel. 

Note that the above argument about the cancellation of the CDW and BCS instabilities 
makes use only of the linearity of the spectrum in a given branch and it is valid for any form 
of the interaction, as long as the around the corner processes can be neglected; moreover, it 
applies also away from half filling. The uncontrolled point of our analysis is the neglect of 
'around the corner' processes, which should be studied more carefully. 

We expect that taking the absence of the around the corner processes for granted, our 
results will stay true in all orders of perturbation theory, since this neglect results in an 
additional conservation law: the number of electrons in a given branch is unchanged in 



the scattering. In D = 1, Dzyaloshinskii and Larkin [[HJ used conservation of electrons 
in a given branch and the corresponding Ward identities to solve the Tomonaga-Luttinger 
model exactly. Recently, the importance of conservation laws in investigating the non-LFL 
behavior in higher dimensions was stressed in several studies [f^|5|,|i"3|| . 

Finally, we solve the model via Haldane's bosonization |Q. In bosonizing the kinetic 
energy H^m = Z) P e pC^ p c p , we closely follow the standard methods in D — 1 (TJ]]. For 
every branch a, we construct a density operator p a (q) = J2 P @«(p + q)©a(p) c P + q c p; where 
B a (p) = 1, if p lies in the quadrant a, and Q (p) = otherwise. Note that D = 2 is different 
from the one-dimensional case in that even restricting ourselves to the neighborhood of the 
Fermi surface, pairs p and p + q lying in different branches exist for small q. This is due to 
the fact that the Fermi surface in D > 1 is a continuous manifold. If we want to interpret 
the physical density operator p(q) as p(q) = J2 a Pa(<i), we have to neglect the excitations of 
the electrons between different branches. This is a justified approximation, if the number of 
neglected processes is much less than their total number, i.e. we have to restrict ourselves 
to the study of sufficiently long wavelength processes, |q| <C A. Neglecting the above 
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mentioned 'around the corner' processes and using the linearity of the spectrum e p , one has 

[H kin , p a (q)] = v Q ■ qp a (q), (8) 

straightforwardly generalizing the result of the one-dimensional calculation. Closely follow- 
ing the derivation of the commutator of the density operators in D = 1 (see, e.g., Ref. 

we find [p Q (q),P/j(-q')] = 6 a ,/3 <W E P (™ p+q - where n° is the distribution of 
noninteracting fermions. From here we have 

[Pa(q),P0(-q')] = -*a,/3 <W 77T77 ™a ' q, (9) 

where n Q is a unit vector in the direction of the group velocity v a . Thus the kinetic energy 
operator can be written H kin = ^ ^ F Xia E q : Pa(q)Pa(~ q) : ; where : X : denotes normal 
ordering of X. Let us investigate the interaction energy, which is given by the second term in 
Eq.(l). Neglecting the around the corner processes, we can write H int = ^ J2 a ,/3 Sq ^(q) : 



p a (q) pp(— q) :• Introducing new operators via p a (q) = a a (q)y ilA/(27r) 2 , the total Hamil- 
tonian can be written 
1 



H = -^2v F \q\ (<W + G) y | cosa cos/3| : a a (q) a^-q) : . (10) 

The commutation relations for a Q (q) are [a a (q), a / g(— q')] = — 5 Q)i g <5 qiq /sign(n Q • q). Thus 
a a (q) is a Bose creation operator for sign(n a • q) = 1 and an annihilation operator oth- 
erwise. We used a, resp. (3 to denote the angles between q and n a , resp. iip. To 
proceed further, we have to choose the direction of q in order to distinguish the cre- 
ation and annihilation operators. It is easy to see that all q, whose angle with the p y 
axis lies in the interval (— 7r/4, vt/4), lead to the same choice of operators. Ignoring the 
ground state energy and concentrating on the excitations, we have after symmetrization 
H' — | J2q V f\ < i\ ^(q)M(q)A(q), where H' is the part of the Hamiltonian corresponding to 
the chosen range of q, A^(q) = (a\, a 3 , a\, a\, a%, «3, a,2, cla) and M is a real symmetric 8x8 
matrix, 



M 



'Mi M 2 ^ 



M 2 Mi 
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The matrices Mi and M2 are given in terms of A = (1 + G)cos0, B = (1 + G)sin0, 



C = G cos (j), D = G sin <p and E = G^J sin cos (where is the angle between ni and q) 
by 

\ 



Mi 



A E 
E B 



M 2 



(ce\ 



E D 







where a 2 is a Pauli matrix and ® stands for a direct product. The Hamiltonian Eq.(lO) 
is quadratic in Bose operators and can be diagonalized by a generalized Bogoliubov trans- 
formation. This problem can be reduced to a standard eigenvalue problem for the matrix 
N = (Mi - M 2 ) (Mi + M 2 ) HI and the Hamiltonian Eq.(lO) can be written 



H = J2v f \ q \ ]T A a &t (q)6a ( q ) 

q q=± 



(11) 



where are the eigenvalues of N, 



1 + 2G ± J(l + AG) cos 2 (20) + (2G) 2 



(12) 



Note that A± as given by Eq.(12) can be defined in the whole q plane by the same formula. 
That is why we omitted in Eq.(lO) the prime on the Hamiltonian and the sum is over the 
whole q plane. In the noninteracting case G = and A± reduce to |sin0| and |cos0|, 
giving the correct excitation energies v a • q for the collective modes. Eqs.(ll,12) represent a 
solution to the model Eq.(l). The low energy spectrum is found to consist of gapless bosonic 
modes. The quasiparticles are dissolved in these collective excitations. 

We would like to point out that it is possible to embed the model spectrum studied 
here into a sequence of spectra e n (k) = (l/2m) (k 2n + fc 2n ) ^ with n — 1,2, ... as its limit 
n — ► 00. In fact, the constant energy line equation e^k 



e has the solution \k r 



V2 



me, 



k y = or \k y \ = v2me, k x — and the spectrum e QO (k) is (in the neighborhood of the 



Fermi energy ep) identical to that studied here, if one requires Vp = ^2ep/m. Decreasing 
n, the Fermi surface becomes more and more curved, until one finally obtains the usual 
isotropic spectrum ej.(k) = k 2 /2m. The crossover between the Luttinger liquid for n = 00 



and Landau's Fermi liquid for n — 1 will be studied elsewhere [17] 



We emphasize that the above analysis does not apply to a half filled nearest neighbor 
tight binding band. In fact, an important point in deriving the cancellation of the CDW and 
BCS instabilities is that the bare spectrum does not depend (within a given branch) on the 
direction perpendicular to v Q , which is not true for a tight binding model. Moreover, in a 
nearest neighbor tight binding band, there are van Hove singularities exactly in the corners 
of the Fermi surface leading to a completely different picture than the one used here. In fact, 
let us consider electrons within the energy shell A around the Fermi surface. In our model, 
the ratio of the number of 'corner' states in the shell to the total number of states within 
the shell is ~ A/vpA —>■ 0, as A — > 0. This justifies the neglect of the corner states in our 
calculation. However, in a nearest-neighbor tight binding model, the fraction of the 'corner' 
states remains finite for A — > and our analysis is not applicable. This is a consequence of 
the van Hove singularities at (±7r,0), (0, ±7r) in that model. Schulz |1| and Dzyaloshinskii 



||T9| applied a complementary method to treat the tight binding spectrum: they showed 
that the BCS and CDW susceptibilities are governed by electrons in the corners of the 
Fermi surface, which led them to replace the Fermi surface by 4 points and they treated this 
simplified model by renormalization group methods. Similarly as in Refs. [T^|l_^] , Houghton 



and Marston |J model a tight binding spectrum close to half filling by 4 points which, 
however, represent the branches of the spectrum and van Hove singularities are neglected. 
In this respect, their model is similar to ours. For spinless fermions with repulsion, Houghton 
and Marston find, in absence of large momentum scatterings, no instability towards CDW 
state formation in agreement with our result. The model Eq.(l) bears some resemblance 



also to the problem studied in Ref. [20], where conclusions similar to ours were reached. 

In conclusion, we studied a two-dimensional model (generalization to higher D is straight- 
forward), where a Luttinger liquid is realized. The model exhibits an analogue of spin-charge 
separation; the role of spin is played by the momentum transverse to the direction of the 
quasi-one-dimensional motion and the electron Green's function is expected to have new 
features compared to D = 1. It will be interesting to study the influence of curvature of the 
Fermi surface, electron spin and large momentum scatterings on the stability of Luttinger 



liquids. 
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FIGURES 

FIG. 1. Fermi surface of the model Eq.(l). The branches are labelled 1-4. A is the length of 
the Fermi surface in one branch. Qi and Q2 are nesting vectors of the Fermi surface. 
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